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Abstract 

We study quenched QCD at finite chemical potential, [ii, for the third component of isospin 
and quenched two-colour QCD at finite chemical potential, fi, for quark number. In contrast to 
the quenched approximation to QCD at finite quark-number chemical potential, the quenched 
approximations to these theories behave similarly to the full theories. The reason is that these 
theories have real positive fermion determinants. In both of these theories there is some critical 
chemical potential above which the charge coupled to the chemical potential is spontaneously 
broken. In each case, the transition appears to be second order. We study the scaling properties 
near the critical point using scaling functions suggested by effective (chiral) Lagrangians and find 
evidence for scaling with mean-field critical exponents in each case. The subtleties associated with 
observing the critical scaling of these theories are discussed. 
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I. INTRODUCTION 



RHIC at Brookhaven and the CERN heavy-ion program give us the possibility of 
producing hadronic matter at high temperatures and finite baryon number density. Cold 
nuclear matter - hadronic matter at finite baryon number density - exists in neutron stars. 
In addition, the bulk properties of large nuclei should be well described as nuclear matter. 
Because of Coulomb interactions, nuclear matter is not only at finite baryon-number density 
but also finite (negative) isospin(/ 3 ) density. Finally, it has been suggested that sufficiently 
dense nuclear matter will also have finite strangeness density. It is therefore of interest to 
study QCD at finite baryon-number, isospin and strangeness densities both at zero and finite 
temperature. 

Such finite densities are achieved by introducing a chemical potential for the rel- 
evant charge operator. Unfortunately, the introduction of a finite chemical potential for 
quark/baryon number leads to a complex fermion determinant which precludes the use 
of standard simulation methods based on importance sampling. Introducing a chemical 
potential fij for the third component (I3) of isospin leaves the determinant real and non- 
negative. Adding the additional (small) I3 breaking term needed to observe spontaneous 
isospin breaking on a finite (lattice) volume, makes this determinant strictly positive and 
simulations possible. Effective Lagrangian analyses of this theory have indicated that it 
should undergo a phase transition at fij = to a state where 1% is spontaneously broken 
by a charged pion condensate with an accompanying Goldstone mode [|I| We are currently 
performing simulations of this theory |3j . Adding, in addition, a finite chemical potential 
(fi s ) for strangeness again makes the fermion determinant complex. In this case, however, 
there are related theories with real positive fermion determinants which mimic the correct 
physics for small /i s 

Until we have a satisfactory way of dealing with finite baryon-number density, it is of 
interest to study models which exhibit some of the anticipated properties of QCD at finite 
quark/baryon-number density but which have real positive fermion determinants. One such 
model is 2-colour QCD at finite quark-number chemical potential (/z). This theory exhibits 
diquark condensation for /i sufficiently large. [|, ^ [| [7], ||, |J Effective Lagrangians and chiral 
perturbation theory analyses of this theory display similar behaviour, and give quantitative 
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predictions for the nature of this transition and the equation of state in the neighbourhood 



of this critical point [jTO], [LI]. These predictions have been validated in the strong coupling 
limit of 2-colour QCD at finite fj, by Aloisio it et al. ||. Formation of such condensates 
has been suggested for QCD at large enough \i JT2], |T3|, [H], [T5], [HJ. There is, of course, 
one crucial difference. The condensate for 2-colour QCD is a colour singlet, the symmetry 
breaking is realized in the Goldstone mode and this theory exhibits superfluidity. For true, 
3-colour QCD, the condensate is, of necessity, coloured, the symmetry breaking is realized 
in the Higgs mode and this theory exhibits colour superconductivity. 

The quenched approximation, i.e. the approximation of setting the fermion deter- 
minant to unity has proved useful for calculating hadron spectra and matrix elements. At 
finite temperatures it shows the deconfinement transition at which chiral symmetry is re- 
stored, as does the full theory. However, since the order of this transition and the equation 
of state are flavour dependent, it yields no useful information on these issues. Where it can 
be used, its principal advantage is that it reduces the computing requirements by several 
orders of magnitude. For QCD at finite quark-number chemical potential, the quenched 
approximation was even more appealing, since it avoids the problem of the complex fermion 
determinant. Unfortunately, it was soon discovered that it does not produce the correct 
physics ]I7|, [L8], [L9|. Whereas it is believed that as fi is increased, the first phase transition 
should occur for fi ~ m N /3, the quenched theory showed a transition for ji rs m n /2. This 
was realized to indicate that the quenched theory should be considered as the zero-flavour 
limit of a theory with an equal number of quark flavours with quark- number +1 and with 
quark- number —1, rather than of a theory where all the quarks have quark number +1. 
This was implicit in this early work |18|, |19| , and was made explicit in terms of random 
matrix models by Stephanov fl2"0| . 

Since QCD at finite isospin chemical potential has equal numbers of quarks with 
I 3 = +| and with J 3 = — |, it thus should be expected to admit a sensible quenched 
approximation. What one immediately realizes is that earlier studies of quenched QCD at 
finite quark-number chemical potential ]I7|, [18], [L9|, [HJ can now be reinterpreted as studies 
of quenched QCD at finite isospin chemical potential. 

Now, however, we can include an explicit I3 violating interaction which makes the 
Dirac operator positive definite (rather than positive semi-definite) giving better conver- 
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gence for our inversion algorithm, and allowing us to measure the /3-breaking pion conden- 
sate directly. Similarly, in 2-colour QCD, since quarks and antiquarks belong to the same 
representation (the fundamental) of SU{2) co i our , it too should have a sensible quenched ap- 
proximation. In fact, we see, it is the same property that gives these theories real positive 
fermion determinants that allows quenched approximations. This is not surprising, since 
all a positive fermion determinant does is to reweight the contributions. This contrasts 
with a determinant where the sign of the real part changes and contributions from different 
configurations can (and in the case of QCD at finite \x must) cancel. 

Once one has determined that each of these theories undergoes a second order transi- 
tion to a state characterized by a condensate which spontaneously breaks the charge coupled 
to the chemical potential, it is useful to examine the scaling properties of the order param- 
eter and certain composite operators in the vicinity of the critical point and to obtain the 
critical exponents. This allows one to write down an equation of state for the system. Such 
equations of state are important for modeling neutron stars. Of course there one needs to 
work at finite baryon-number chemical potential as well. Effective Lagrangians and chiral 
perturbation theory through 1-loop suggest that we should see mean field scaling near this 
critical point % |TD|, [□]. 

We have measured the pion condensates, chiral condensates and isospin densities as 
functions of isospin chemical potential (/xj) and the explicit isospin breaking parameter A on 
12 3 x 24 and 16 4 quenched QCD gauge configurations at (3 = Q/g 2 = 5.7, and on 8 4 quenched 
gauge configurations at (3 = 5.5. Both show evidence for the expected mean-field scaling. 
In addition we have measured the diquark condensate, chiral condensate and quark-number 
density on a set of 8 4 and 12 4 quenched 2-colour QCD configurations with (3 = 4/g 2 = 2.0. 
Here mean-field scaling is again favored. 

In section 2 we present our actions and their relevant symmetries. Section 3 describes 
our simulations and results. Critical scaling analyses are presented in section 4. Discussions 
and conclusions are presented in section 5. 
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II. ACTIONS AND SYMMETRIES 



The staggered fermion part of the action for lattice QCD with degenerate u and d 
quarks at a finite chemical potential /xj for isospin (J 3 ) is 



(1) 



sites 



where p(n) is the standard staggered p with links in the +t direction multiplied by 
and those in the —t direction multiplied by e~5^. The explicit symmetry breaking term 
i\e~x r T2X) the lattice equivalent of iXip'j^ip, is in a direction in which the J3 symmetry is 
expected to break spontaneously for fij sufficiently large. This term is necessary in order 
to observe spontaneous symmetry breaking from a finite lattice. We will be interested in 
the limit A — > 0. We will present a detailed discussion of the symmetries of this theory in a 
forthcoming paper on our simulations with dynamical quarks. The Dirac operator 



M 



Pi/J-i) +m Ae 

-Ae p(-fi I )+m 



(2) 



has determinant 

det M = det [A*A + A 2 ] (3) 

where we have defined 

A=p{m)+m. (4) 

We note that this determinant is positive for A 7^ 0, as promised. Observables we measure 
include the chiral condensate 

(W) & ( xx ), (5) 

the charged pion condensate 

i(4>lbT 2 ^) & i(xtT 2 x) 

and the isospin density 



1 dS 



f 



The quark action for 2-colour QCD with one staggered quark is 
Sf = E \x[P(») + m] X + \^\X T T2X + XT2X T ]\ 

sites 1 A ' 



(6) 
(7) 



(8) 
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where is the normal staggered covariant finite difference operator with \x introduced 

by multiplying the links in the +t direction by e M and those in the —t direction by e" 1 '. The 
superscript T stands for transposition. The term proportional to A explicitly breaks quark- 
number symmetry, and again we shall be interested in the limit A —>■ 0. The symmetries of 
this action and the positivity of the determinant and pfaffian have been discussed extensively 
in previous work on dynamical quark simulations and will not be repeated here 0, |5], |7J . 

The observables we measure include the chiral condensate, 

(xx) = (W>), (9) 

the diquark condensate, 

(X T T2X) (10) 

and the quark-number density 

III. QUENCHED SIMULATIONS AT FINITE CHEMICAL POTENTIALS 
A. Quenched QCD at finite isospin chemical potential 

We have calculated quark propagators at a finite chemical potential \ii for the third 
component, J3, of isospin from a noisy source on an ensemble of 12 3 x 24 and 16 4 quenched 
gauge field configurations at f3 = 5.7 and on 8 4 quenched gauge field configurations at 
(3 = 5.5. From these propagators we obtained stochastic estimators for the pion condensate, 
the chiral condensate and the isospin density as functions of fij and the ^3 breaking parameter 
A. We chose A << m, since we are interested in the limit A — > 0. 

Let us first consider the (3 = 5.7 calculations. We generated 100 (3 = 5.7 
12 3 x 24 equilibrated configurations separated by 1000 sweeps consisting of 9 overre- 
laxation sweeps followed by 1 10-hit metropolis sweep repeated 100 times, which ap- 
peared to generate relatively independent configurations. On each of these configura- 
tions we obtained a stochastic estimate of i(X eT 2X)i (Xx) an d Jo using a single noisy 
source, for m = 0.025, A = 0.0025,0.005,0.0075 (and A = for m < /i c ) and 
Hi = 0.0,0.1,0.2,0.3,0.35,0.4,0.45,0.5,0.55,0.6,0.7,0.8,0.9,1.0,1.2,1.6,1.8,2.0,2.1. The 
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fact that we used the same set of configurations and noise vectors for each value of A and hi 
means that all our 'data' points are strongly correlated. We also measured the pion mass at 
/i = A = using a wall source, on these configurations gauge fixed to Coulomb gauge. We 
obtained m n = 0.441(1). Runs at A = 0.0025 were performed on an 8 4 lattice at selected Hi 
values to check that finite size effects were not too large. 

From these 12 3 x 24 measurements we determined that the neighbourhood of 
Hi = 0.45 was of particular interest. We then increased our spatial volume and per- 
formed simulations on a 16 4 lattice, again at f3 — 5.7 and m = 0.025, to allow us to 
run at smaller A. Here we generated our quenched configurations and performed our mea- 
surements of the condensates and isospin density using the same hybrid molecular-dynamics 
code that we use for simulations with dynamical quarks. This meant that we used an 
independent set of configurations for each value of A and Hi- 100 independent config- 
urations separated by 10 molecular-dynamics time units were generated for each A and 
Hi- (With no fermions in the updating, 10 time units were adequate to decorrelate con- 
figurations. The time increment for updating was dt = 0.1 which is adequate without 
light fermions in the updating.) We chose A = 0.000625,0.00125,0.001875,0.0025, and 
used fii = 0-2, 0.25, 0.3, 0.35, 0.375, 0.4, 0.425, 0.45, 0.475, 0.5, 0.525, 0.55, 0.6, 0.65, 0.7 to ad- 
equately cover the region of interest. 

The charged pion condensates are shown in figure |l| What is immediately clear is that 
the condensate is small for small Hi- At /i/ = 0, the condensate vanishes as A — > 0. Aloisio 
et aZ.|| have shown that this vanishing of the condensate extends to small but finite hi- 
We have performed various extrapolations to A = — linear, linear plus cubic, linear plus 
quadratic - all of which are based on the predictions of partially conserved current analyses 
and their implementations in effective (chiral) Lagrangian models. These extrapolations 
suggest that this condensate vanishes out to at least hi ~ 0.4. (the A = measurements are 
trivially zero and do not enter into this discussion). More importantly, these extrapolation 
schemes, and those based on power law scaling in A as would be expected for hi at a critical 
value, indicate that it is highly probable that this condensate remains finite as A — > in the 
infinite volume limit for \xi ~ 0.5. This indicates that there is a phase transition for A = at 
some fii — He with 0.4 < hi ~ 0.5, below which the pion condensate is zero and above which 
there is a non-zero charged pion condensate which breaks J3 spontaneously with associated 



7 



SU(3) quenched 0=5.7 m=0.025 12 3 x24 lattice 
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SU(3) quenched = 5.7 m=0.025 16 4 lattice 
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FIG. 1: a) Charged pion condensate, i{ip r Y5T2ip) = z(x er 2X)> as a function of fij, on a 12 3 x 24 
lattice with [5 = 5.7 and quark mass m = 0.025. b) Charged pion condensate on a 16 4 lattice. The 
solid lines superimposed on our measurements are the scaling fits described in section 4. 
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Goldstone pions. Since there is no abrupt jump at /// = /i c , this suggests that the phase 
transition is second order. This condensate increases monotonically up to \ij <^ 1.6 after 
which it falls rather rapidly approaching zero by fij ~ 2.0. This fall is clearly a saturation 
effect, a conclusion that becomes even more compelling when we look at j'q. 

SU(3) quenched £ = 5.7 m=0.025 
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FIG. 2: The chiral condensate, (ipip) = (xx)-> at /3 = 5.7 and m = 0.025 as a function of [ij. 

Figure shows the chiral condensate as a function of \ii for both the 12 3 x 24 and 
16 4 simulations at (3 — 5.7, m = 0.025. We see that it remains approximately constant at 
its fi = values for fij < fi c , above which it commences its descent towards zero. What is 
also clear when looking at the plots of chiral and pion condensates, is that the expectation 
from tree level chiral perturbation theory that the condensate merely rotates from the chiral 
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symmetry breaking to the isospin breaking direction has at best a rather limited range of 
validity. 
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FIG. 3: Isospin density at /? = 5.7 and m = 0.025 as a function of /i/. The curves are from the 
scaling fits to the pion condensate, and are described in section 4 

The isospin density is presented in figure ^. The main features are its slow rise 
from zero above /x c , followed by a more rapid rise at larger fij, finally flattening out at its 
saturation value of 3 for \ij > 2.0. Note that the saturation value is 3 because we have 
chosen to normalize it to 8 continuum flavours, the number associated with the action given 
in section 2. The value 3 arises from | for each of 3 colours and 2 staggered 'flavours' (= 8 
continuum flavours) on each site, the maximum allowed by fermi statistics. (For contrast 
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the condensates have been normalized to 4 flavours for comparison with those reported in 
finite temperature simulations.) 

We now turn to the consideration of our simulations and measurements at stronger 
coupling, (3 = 5.5, which were performed on an 8 4 lattice. Again we chose m = 0.025 in 
lattice units. Here we followed the same procedure as for the /3 — 5.7, 16 4 simulations, using 
the same hybrid molecular-dynamics code as used for our dynamical simulations. We used 
A = 0.0025, 0.005, 0.0075 (and A = for fij < /i c )in lattice units. Our chosen /// (in lattice 
units) were 0.0, 0.1, 0.2, 0.3, 0.35, 0.4, 0.45, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.2, 1.4, 1.6, 1.8, 2.0, 2.2. 
We checked for finite size effects, by repeating these simulations on a 12 4 lattice for se- 
lected \ij values (0.3, 0.4, 0.5, 0.6, 0.8 and 1.0), with A = 0.0025 where such effects are 
expected to be largest, again analyzing 100 configurations for each /i/. For each /i/ the 
difference between the 8 4 and 12 4 measurements was within 2 standard deviations of zero. 

The /3-breaking charged pion condensate is plotted in figure f| as a function of /i/, for 
the A values mentioned above. Superficially, this graph appears similar to figure [L[ Again, 
we note that any reasonable extrapolation to A = would suggest that, in this limit, the 
condensate vanishes for fij < 0.3, while it clearly does not for /i/ > 0.5, which strongly 
suggests that there is a phase transition at /i/ = /i c with /i c lying in the range 0.3-0.5. (We 
have not performed such extrapolations since at larger fij values, statistical fluctuations lead 
to incorrect ordering of the condensates for different As, which would produce unphysical 
extrapolations.) However, what one notices on closer inspection is that the transition is 
somewhat steeper than that for (3 = 5.7. We will return to this point in the next section. 
Saturation again sets in for \xi a little above 2.0. 

The chiral condensate shown in figure [| behaves very similarly to that at (3 = 5.7. 
The main reason for presenting it here is to display the scaling predictions of section 4. We 
note, however, that its value is 0.764(7) at A = /ij = 0, so again it is at most over a limited 
range of /zj that the condensate simply rotates from the chiral towards the isospin-breaking 
direction. However, the violations are less severe than at (3 = 5.7. The isospin density 
also behaves very similarly to that at weaker couplings, rising steadily from zero as \x\ is 
increased past /z c and eventually saturating at 3. This is plotted in figure || again, primarily 
to show the scaling predictions described in the next section. 
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SU(3) quenched 0=5.5 m=0.025 8 4 lattice 
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FIG. 4: Charged pion condensate as a function of fij and A on an 8 4 lattice with j3 = 5.5 and 
m = 0.025. The solid lines are fits described in section 4. 

B. 2-colour QCD at finite quark-number chemical potential. 

We calculated the quark propagators for fundamental staggered quarks at fi- 
nite chemical potential fi from a single noisy source for each of a set of quenched 
SU(2) gauge field configurations at (3 = 2.0. These configurations were generated 
using the hybrid molecular-dynamics code of our dynamical quark simulations. We 
generated 100 such independent gauge configurations spaced by 10 molecular-dynamics 
time units for each \x and A value. Our simulations used dt = 0.1 for the up- 
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SU(3) quenched 0=5.5 m=0.025 8 4 lattice 
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FIG. 5: The chiral condensate as a function of /// for quenched QCD at j3 = 5.5, m = 0.025. The 
solid lines are fits derived from the scaling analyses of section 4. 

dating. We used A = 0.0025,0.005,0.0075,0.01 (and zero for fi < fi c ), for // = 
0.0, 0.1, 0.15, 0.175, 0.2, 0.21, 0.225, 0.25, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.95, 1.0, 1.1, which ad- 
equately covered the range of interest. Since the 8 4 'data' showed evidence of finite size 
effects, we repeated these simulations on a 12 lattice. From these noisy quark propagators 
we obtained stochastic estimators for the diquark condensate, the chiral condensate and the 
quark- number density, each normalized to 4 continuum flavours (1 staggered quark field). 

Figure |7| shows the diquark condensate as a function of fi for each of the considered 
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SU(3) quenched 0=5.5 m=0.025 8 4 lattice 
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FIG. 6: The isospin density as a function of [ii for quenched QCD at f3 = 5.5, m = 0.025. The 
solid lines are fits derived from the scaling analyses of section 4. 

A values. For fi < 0.15 it is clear that any reasonable extrapolation to A = from finite 
A would yield estimates close to zero for the condensate in this limit. For /x > 0.25, the 
extrapolated condensate clearly remains greater than zero for A — > 0. Thus we deduce 
that there is a phase transition at /j = /j c from the normal state to one in which quark 
number is broken spontaneously by a diquark condensate, for some \i c between 0.15 and 
0.25. Such a symmetry breaking would be accompanied by a diquark Goldstone boson. 
The reason we have not performed an extrapolation in A is clear when one looks at the 
behaviour of the condensate at A = 0.0025 relative to the higher lambdas. We note that 
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SU(2) quenched ,8 = 2.0 m = 0.025 8 4 lattice 




SU(2) quenched ,8 = 2.0 m=0.025 12 4 lattice 
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FIG. 7: The diquark condensate for quenched 2-colour QCD as a function of fi and A on at (3 = 2.0 
a) on an 8 lattice and b) on a 12 4 lattice. The curves are the scaling fits described in section 4. 

this condensate increases monotonically with /i up to fi ~ 0.7 after which it falls due to 
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saturation effects. Qualitatively this curve is very similar to the corresponding curve with 
dynamical quarks. It is also very similar to the curve for QCD at finite isospin chemical 
potential at (5 = 5.5 especially when one observes that we should equate fii with 2/i. As in 
that case, the condensate increases relatively rapidly as [i is increased beyond [i c . (This will 
be made quantitative in the next section.) We note that the A = 0.0025 measurements on 
the 8 4 lattice show signs of finite size effects for larger /i values, which is why we repeated 
these simulations on a 12 4 lattice where these effects have been considerably reduced. 
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SU(2) quenched 0=2.0 m=0.025 12 4 lattice 
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FIG. 8: Chiral condensate as a function of \x and A, for 2-colour QCD at finite chemical potential 
for quark number on a 12 4 lattice. The solid lines are the fits given in the next section. 

In figure |8| we show the chiral condensates as functions of /i for the 4 A values given 
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above (and A = for /x < 0.2). (We present only the 12 4 'data' since there is little difference 
between the 2 lattice sizes.) For /x < /x c this graph is consistent with the expectation that 
the A — > limit should be /x independent. As \i is increased beyond /x c , the condensate 
falls. This agrees with the expectation that the condensate is rotating from the chiral to the 
diquark direction in this region. One notes, however, that this is not a simple rotation since 
the magnitude of the condensate also increases with /x for larger /x values until saturation 
effects appear. The general appearance of this curve is similar to those with dynamical 
quarks and to those for QCD at finite isospin density. 

SU(2) quenched 0=2.0 m=0.025 12 4 lattice 
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FIG. 9: Quark-number density jq as a function of zx and A on 
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The quark-number density is shown in figure ^ as a function of /i for each of the 
considered A values mentioned above. (Again, because the finite size effects are small, we 
present only the 12 4 'data'). We see that jo is consistent with zero for small /i, and rises 
slowly from zero above /i c . The rise becomes more rapid for fi > 0.6. The density finally 
saturates at 2 — 2 colours of staggered fermions per site — due to fermi statistics, a finite 
lattice spacing effect. Note that there is little if any A dependence. This behaviour is very 
similar to that of j'q for our quenched QCD at finite /xj, except that there the saturation value 
is 3 rather than 2 since this number merely counts colours. Again we note the similarity of 
these quenched results to those obtained with dynamical quarks. 

The qualitative behaviour of all the 'data' presented in this section is in agreement 
with what is predicted from effective Lagrangians using the methods of chiral perturbation 
theory. We have noted that the condensate does not simply rotate from the chiral direction 
into the quark-number or isospin violating direction which was the prediction of tree level 
chiral perturbation theory [1], lQfl , but rather rotates and increases in magnitude. Such 
behaviour is seen when the chiral perturbation theory includes next-to-leading order terms 
including 1-loop contributions pj| . Of course, the one effect that effective Lagrangians 
cannot see is saturation, since they do not describe the underlying fermions of the theory. 
However, as we have pointed out, saturation is a finite lattice spacing effect, and is thus of 
limited interest. 



IV. CRITICAL SCALING AND THE EQUATION OF STATE. 

Analyses of effective Lagrangians using the methods of chiral perturbation theory 
give predictions for the critical exponents and the scaling of the order parameter(s) close to 
the critical point /i c , both for 2-colour QCD at finite fi and for QCD at finite [ij. Clearly 
the tree level analyses will give mean-field scaling with critical exponents (3 m = ~ and 5 = 3 
]T|, |TD[ . More recent analyses have extended these results to next-to- leading order including 
1-loop corrections for 2-colour QCD, and find that the critical indices remain at their mean- 
field values [JO]]. However, the tree- level result that the magnitude of the total condensate is 
independent of /i and A no longer holds. One can thus speculate that scaling with mean-field 
exponents holds to all orders. The effective Lagrangian for QCD at finite fii is so similar 
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as to lead to the expectation that this theory also has mean field behaviour through next- 
to-leading order and perhaps to all orders. Simulations of 2-colour QCD at finite \x in the 
strong coupling limit |8j also indicate mean-field scaling with a scaling function consistent 
with that of lowest order tree-level effective Lagrangian predictions. 

The two critical exponents relevant to the scaling of the charged pion condensate - 
the order parameter for the transition we are considering — are /3 m and 8 defined as follows. 
At m = 0, 

i{i>lsT 2 V) -> const(m - n c f m (12) 
as hi — > fM c + (and vanishes for fii < /-O- For fii = A^c 

i(' i l ! l5 T 2'4 ! ) — > const m* (13) 

as m — > 0+. The mean- field values of these 2 exponents are (3 m = \ and 8 = 3. Analogous 
definitions hold for 2-colour QCD. 

We shall be considering 3 different forms for the scaling function, i.e. for the equation 
of state, each of which exhibits mean-field scaling. Indeed, sufficiently close to the critical 
point they give identical scaling functions for the order parameter. Our later discussions 
will make it clear why, in practical applications, it is important to have a more detailed 
knowledge of the scaling function. 

The first form for the equation of state is that suggested by the lowest order tree-level 
analysis of the effective Lagrangians of JjJ [ID| , which are of the non-linear sigma model form. 
It derives from the value of the parameter a which minimizes the effective potential 

£ = —a /i 2 sin 2 (a) — bm cos(a) — b A sin (a) (14) 

in terms of which 

(X T T2X) = b sin(a) (15) 
(XX) = b cos(a) (16) 

and 

jo = 2 a [A sin 2 (a). (17) 
b is given in terms of /x c and a, namely 

2 

b = — au 2 r (18) 
m 
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and the scaling forms of the observables for A = 0, expressed in terms of the scaling variable 
x — nl /i c , are 

(X T T2 X ) = b ^1-^ (19) 
(XX) = b ^ (20) 

and 

j = 2 OA *(l-i) (21) 

Note that the equivalent form for QCD at finite /i/ has an extra factor of 2 in the equation 
for j'q because of our normalization. Otherwise the forms are identical with fi replaced 
with fij. The main deficiency with this form is that it requires that the magnitude of the 



condensate, J (x Tr 2X} 2 + {xx} 2 to be constant independent of /x and A. This is not even 



true in next-to-leading order chiral perturbation theory [] L] and has at best a limited range 
of validity in our simulations. 

To allow the norm of the condensate to increase with // and possibly A, we consider 
a second form for the equation of state which is derived from an effective Lagrangian of 
the linear sigma model type. This is obtained by extracting the values of R and a which 
minimize the effective potential 

11 1 

£ — - R 4 — - a R 2 — —b fi 2 sin 2 (a) — cm R cos(a) — cXR sin(a) (22) 



in terms of which 



and 



c is given in terms of \x c by 



(X T T2X) = c R sin(a) (23) 

(XX) = cR cos(a) (24) 

j = b^R 2 sm 2 (a). (25) 
b\i\ 



m 



a + bjj? c (26) 



When A = we obtain the scaling forms 



(x T T2X) = c ] ja + b ^ - (a + b fil) 1 (27) 
(Xx) = cyja + bti 2 ± (28) 

Jb 
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and 

jo = b fi a + b fi^x 2 - (a + b fjfy—^ . (29) 
Again j'q has an extra factor of 2 for QCD at finite ///. 

Finally we consider the simplest mean-field equation of state |22| which considers only 
the scaling of the order parameter (diquark or charged pion condensate) near the transition 
and does not even attempt to address chiral symmetry breaking. This is given in terms of 
cf) which minimizes the effective potential 

8 = -\ a (v-Vc)<l> 2 -b\<t>. (30) 

In terms of <f), 

{x T T2X) = b(t> (31) 

and 

Jo = a 2 . (32) 

Once again Jq has an extra factor of 2 for QCD at finite ///. When A = we find the simple 
scaling relations for /x > \i c 

(X T T2X) = byja(/j, - fi c ) (33) 

and 

jo = a 2 (/j, - /i c ). (34) 

We start by considering quenched QCD at finite /// at (3 = 5.5. We first fit our 
measurements of the pion condensate to the form suggested by the non-linear sigma-model 



effective Lagrangian (equations ^,15). The quark mass is treated as a fitting parameter. 
Here we find a good fit over the range < /// < 0.7, and all 3 A values with \i c = 0.422(2), 
a = 0.0618(8) and m = 0.0267(2) with a confidence level of 61%. The range limitation 
is due to the fact that the norm of the condensate is only approximately constant over a 
limited range of /i/. The departure of the mass m from its input value is presumably also a 
reflection of this. We then fit our pion condensate 'data' to the form suggested by the linear 
sigma-model effective Lagrangian (equations ^,^). With this form, our fits constrain the 
mass parameter, m, to be so close to the input mass that we choose to set m = 0.025 in our 
fits. We then obtain an excellent fit over to the pion condensate over the range < /// < 1.4 
with fi e = 0.419(2), a = 0.412(7), b = 0.242(4) and a confidence level of 92%. We have 

21 



superimposed these fits on the graphs of our 'data' in figure |[ These fits give predictions 
for the chiral condensate and isospin density (equations p4T , p5|) which we have included in the 
plots of our measurements of these quantities in figures |||| The predicted chiral condensate 
is clearly in good agreement with our measurements. The predicted isospin density is in 
good agreement with measurements up to fij ~ 1.0. Note that the isospin density has been 
predicted to be proportional to /xf at large /i/ |T| so that the predictions of equation ^5] are 
guaranteed to break down, even if we did not have the effects of saturation. 

Let us now turn our attention to our simulations of 2-colour QCD at finite /i. First 
let us consider fits to our 8 4 measurements. Here we again note that our fits to the non- 
linear sigma model scaling form can only be made over over a somewhat limited domain of 
/i. This can be extended by using the linear sigma model scaling forms, but only when we 
neglect the A = 0.0025 measurements which show clear signs of finite size effects at large 
/i. Fitting the A = 0.005,0.0075,0.01 'data' we obtained a fit over the range < /i < 0.7 
with /i c = 0.209(1), a = 0.405(9), b = 0.71(2) and a 77% confidence level. To reduce the 
finite size errors, especially at A = 0.0025, we repeated these simulations on a 12 4 lattice. 
Here our best fit, obtained using all 4 A values and fitting over the range < ji < 0.4, gave 
/i c = 0.2066(4), a = 0.405(9) and b = 0.73(1) at a confidence level of 10% (x 2 /dof = 1.3). 
Although we have reduced our finite size effects by going to a larger lattice, we have also 
reduced the statistical errors, making the remaining finite size effects more important. We 
suspect that this is at least part of the reason that we have obtained a worse fit on the larger 
lattice. These fits for the pion condensate are shown in figure |7|. The predictions for the 
chiral condensate and quark-number density which these fits provide are superimposed on 
the plots of measured values for the 12 4 lattice in figures [8]||. It is clear that the prediction 
for the chiral condensate provides good fits to the 'data'. The predicted quark-number 
density is also a good fit to the 'data' for fij < 0.4. 

Finally we turn to the consideration of the (3 = 5.7 quenched QCD simulations at 
finite fij. After producing 'data' over the range < \ij < 2.1 for A = 0.0025,0.005,0.0075 
on a 12 3 x 24 lattice, we performed simulations on a 16 4 lattice over the range 0.2 < 
/// < 0.7 at A = 0.000625,0.00125,0.001875,0.0025 to probe the transition more closely. 
Not surprisingly, considering that the norm of the condensate grows rapidly in the broken 
symmetry phase, fits to the non-linear sigma model scaling proved fruitless. We fit our 
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16 4 measurements the linear sigma model equation of state, treating the quark mass as a 
fitting parameter. We exclude the A = 0.000625 measurements from our fits since they 
show clear finite size effects at larger \xj. Our 'best' fit to the A = 0.00125,0.001875,0.0025 
measurements gives li c = 0.4500(7), a = 7.00(1), b = 0.266(7), m = 0.0212(4) with x 2 /dof = 
1.7 and a confidence level of only 0.3%. This indicates that while the linear sigma model 
equation of state remains a good guide to the critical behaviour quenched QCD at finite ixj at 
(3 = 5.7, the data is beginning to show small, but non-negligible departures from this simple 
form. We have superimposed these fits on the 'data' of figure [l] to indicate that they are a 
good qualitative fit to the data. We note also that the predicted fx c is close to the measured 
pion mass m n = 0.441(1). Because m ^ 0.025 in these fits, we do not expect the predictions 
for the chiral condensate to be good. In fact, what we see is that the fits are suppressed 
relative to the data roughly in the ratio of the fitted mass to the input mass (0.025). Next we 
plot the values of jg predicted using equation ^5] to compare with measurements in figure |3|. 
We see that there is good agreement out to \xi = 1.0. Finally, we fit the pion condensate at 



ix i = 0.45 which is very close to fx c to the scaling form of equation [L3|. For this we use the 
measured values on the 16 4 lattice at A = 0.000625, 0.00125, 0.001875, 0.0025 and those from 
the 12 3 x 24 lattice at A = 0.005, 0.0075. From this we obtain 5 = 3.25(7) with a confidence 
level of 27%. Dropping the 2 largest A values gives 5 = 2.88(16) with an 85% confidence 
level. These fits are in good agreement with the mean field value 5 = 3. 

The effective Lagrangian approach seeks to give a uniform description of both chiral 
symmetry breaking and isospin (or quark-number) breaking. Thus we should only expect 
physics to be insensitive to the details of the effective Lagrangian/effective potential when 
the quark mass is small enough that the order m or equivalently order m\ corrections are 
small. One measure of this is how much the chiral condensate (ipip) at the quark mass we 
use differs from its zero mass value. At (3 = 5.5, {ipip){m = 0.0125) = 0.729(11), (ijjtfj)(m = 
0.025) = 0.764(7) and (^)(m = 0.05) = 0.838(5), so the difference is ~ 10% which should 
be small enough. At = 5.7, (#)(m = 0.00625) = 0.288(3), (#)(m = 0.0125) = 0.323(2), 
(ifjif))(m = 0.025) = 0.386(2) and (ifjip)(m = 0.05) = 0.494(2), so the difference is ~ 50% 
and the detailed form of the effective potential might be expected to be important. We 
believe that this is for this reason why the simple ansatz for the effective potential/equation 
of state works so well at /3 = 5.5 but shows limitations at (3 = 5.7. 
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We end this section with a discussion of the behaviour of different variables which 
describe the scaling in terms of For convenience we express them in terms of x = 

H/ fi c (fii/ n c ). The simplest is x — 1 which is the scaling variable of the simplest mean-field 
equation of state. The next, |(1 — 1/x 4 ), is that of the non-linear sigma model (or the linear 
sigma model in the limit a — > oo). The third, |(x 2 — 1/x 4 ) is that for the linear sigma model 
in the limit as a — > 0. This third form gives a reasonable description of the scaling of (5 — 5.7 
'data' on the 12 3 x 24 lattice, in the ordered (fij > /i c ) phase. All 3 scaling variables behave 
identically as x — > 1, however, a poor choice leads to a severely curtailed scaling window. 
x - 1 and \{l - 1/x 4 ) already differ by ~ 10% at x = 1.04. x - 1 and \(x 2 - 1/x 4 ) differ 
by ~ 10% at x = 1.08. Thus, for practical reasons, we should make a wise choice of scaling 
variables. 

However, we have noted apparent mean-field scaling of our (5 = 5.7 pion condensate 
in the ordered phase both in terms of the simplest mean-field scaling forms where the scaling 
variable is x — 1, and in terms of the linear sigma model form where the scaling variable is 
approximately |(x 2 — 1/x 4 ), out to x ~ 1.8. Only by examining the scaling very close to 
the transition and including /// < \x c measurements in the fit, were we able to determine 
that the linear sigma model form gave the better description of the transition. The reason 
for this is that \(x 2 — 1/x 4 ) has a point of inflection at x — \/l0 = 1.467..., at which it is 
again linear in x. In the neighbourhood of this point of inflection |(x 2 — 1/x 4 ) oc x — x , 
where x = (5/8)v^l0 = 0.917..., so this defines another larger region where ^(x 2 — 1/x 4 ) is 
linear in x, and the linear approximation has an intercept not too far from the true critical 
point. It is this linearity, which gives us an (apparent) extended scaling regime. Over the 
region where the scaling variable is approximately linear in x, we will observe scaling with 
the correct exponent, f3 m . Of course, the coefficient of this apparent scaling form will differ 
from the correct coefficient. 

It is instructive to also look at the properties of |(1 — 1/x 4 ) as a function of x. It 
clearly does not have a point of inflection, however, its square has a point of inflection at 
x = ^9/5 = 1.158... about which it too becomes linear. This linear approximation has its 
intercept at x = (9/10) -^9/5 = 1.042..., which is again near to the critical point. Since it 
is the square of the scaling variable that has this point of inflection, this leads to apparent 
scaling with a critical exponent, f3 m of half the true value. In the case of a mean-field critical 



24 



point this pseudo-exponent will be ~, which is [3 m for a tricritical point. If we restrict 
ourselves to the ordered phase, it is possible to fit the order parameter for QCD at finite Hi 
at j3 = 5.5 and that for 2-colour QCD at finite \i at /3 = 2.0 to a tricritical scaling form - 
equation |30] with replaced by ■ Hence one should be suspicious if an order parameter 
scales with a critical exponent which is half of what one expects. 



V. DISCUSSION AND CONCLUSIONS 

Gauge theories with fermions at finite chemical potential for a conserved charge, 
which have positive fermion determinants, have sensible quenched approximations. This 
contrasts with QCD at a finite chemical potential for baryon/quark number, where the 
quenched limit is really the quenched limit for QCD with equal numbers of quarks with 
quark-number ±1 and not that for QCD with all quarks having quark-number +1. 

We have studied the quenched versions of QCD at finite chemical potential (/ij) for 
isospin, and of 2-colour QCD at finite chemical potential (/i) for quark number. We included 
an explicit symmetry breaking term with strength A in each case. What we find is that for 
A — > this symmetry — I 3 for QCD at finite fij, quark-number for 2-colour QCD — is 
unbroken for small fi(fii). In each case there is a critical value for fij or fi (fi c say), above 
which a condensate forms which breaks the relevant symmetry spontaneously. In our QCD 
at finite /i/ measurements at (3 — 5.7 on a 12 3 x 24 lattice, we also measured the pion mass 
at fij = A = 0, and found it to be close to our estimates for /i c as expected. In QCD at 
finite \ii this is a charged pion condensate while for 2-colour QCD it is a diquark condensate. 
This is precisely the behaviour expected and seen for the full (unquenched) theories. The 
transition appears to be second order as was observed and expected in the full theories. For 
high fii(fi) saturation effects due to fermi statistics drive the condensates to zero which is 
also observed for the full theories and is a finite lattice spacing effect. Thus these quenched 
theories are useful to study since they require far less computing time than the full theories 
- these studies were performed on 400MHz PC's supplemented with a few days running on 
an IBM SP whereas the full theories require many months of IBM SP running. 

We studied the scaling of the order parameter, i^'y^ip) (the charged pion conden- 
sate), for QCD at finite fii, and (x Tt 2X) (the diquark condensate) for 2-colour QCD at 
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finite /i and saw evidence that the mean field equation of state was the correct one in the 
neighbourhood of the critical point. This required using scaling forms implied by effective 
Lagrangian analyses, in order to use scaling variables which maximize the scaling window 
and make scaling analyses practical. Because of this sensitivity to the choice of scaling vari- 
ables, these conclusions cannot be considered conclusive. In the quenched theory, it might 
be feasible to extend these simulations to much larger lattices and much smaller A values, 
enabling one to probe sufficiently close to the critical point to be insensitive to the choice 
of scaling variables. However, here one runs the risk of being frustrated by the pathologies 
of quenching. For the dynamical theories the expense (in computing resources) of simu- 
lating on large lattices with very small A values is prohibitive. We have noted that even 
by couplings as small as that at f3 = 5.7, the simple scaling forms suggested by effective 
Lagrangian analyses are barely adequate at the quark mass we use (m = 0.025 in lattice 
units), and a better choice is desirable and probably essential at weaker couplings. While it 
is easy enough to produce variant forms of the effective potential and the equations of state 
they imply, one needs a form which improves the fits, without introducing too many extra 
parameters and preferably one motivated by physics or chiral perturbation theory. 

On examining the detailed properties of proposed scaling variables, we have found 
that a certain class has a point of inflection, when expressed in terms of the simplest scaling 
variable. This leads to an extended region of apparent scaling which is a precursor of true 
scaling. In dynamical theories where one cannot probe the true scaling window, this might 
be the best indicator we have of the nature of critical scaling. However, as we have noted, 
certain other proposed scaling variables which lack such a point of inflection, have squares 
which do have such a point of inflection, leading to apparent scaling with a critical index 
f3 m which is half of the true value, potentially leading one to erroneous conclusions. 

Once the fit to the pion or diquark condensate provides the parameters for the appro- 
priate effective potential, this provides a prediction for the isospin or quark-number density. 
These predictions are in good agreement with our direct measurements of these 'charge' den- 
sities within the scaling window. In addition, the effective potential provides a prediction 
of the chiral condensate. For quenched QCD at finite /// at (3 — 5.5 and quenched 2-colour 
QCD at finite fj, at (5 — 2.0, where the fits require a quark mass m in good agreement with 
the input value (m = 0.025), the predicted chiral condensates are in good agreement with 
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the measured values. In quenched QCD at finite fij at (3 = 5.7, where the fits require m to be 
appreciably lower than the input value, the predictions for the corresponding chiral conden- 
sate are similarly lower. This we tentatively attribute to the fact that the higher order terms 
in a chiral expansion (expansion in powers of m) are large at f3 — 5.7, m = 0.025. Where the 
predictions from our proposed equation of state are in good agreement with measurements, 
they support our claim that there is a critical point at = fi c , A = 0, with mean field 

exponents. 

We have just completed further simulations of 2-colour QCD at finite quark-number 
chemical potential with dynamical quarks which examine the critical scaling at the transition 
from the normal phase to the diquark-condensed phase in which we also see evidence for 
mean-field scaling [p6|| . In addition, we are performing simulations of QCD at finite chemical 
potential for isospin(/ 3 ) |2], |3|] and see evidence for mean field scaling 0. 

We are extending the work of this paper to include a strange quark with its own 
chemical potential. This work will study the competition between charged pion and kaon 
condensation as the 2 chemical potentials are varied independently. Chiral perturbation 
theory predicts that the phases with pion and kaon condensates are separated by a first 
order phase transition. Here the quenched approximation will make it easy to study the 
phase structure of the plane defined by the 2 chemical potentials. 
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